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Abstract 

In this paper we prove the existence of vortices, namely standing waves with non null angular 
momentum, for the nonlinear Klein-Gordon equation in dimension N > 3. We show with vari- 
ational methods that the existence of these kind of solutions, that we have called hylomorphic 
vortices, depends on a suitable energy-charge ratio. Our variational approach turns out to be 
useful for numerical investigations as well. In particular, some results in dimension N = 2 are 
reported, namely exemplificative vortex profiles by varying charge and angular momentum, to- 
gether with relevant trends for vortex frequency and energy-charge ratio. The stability problem 
for hylomorphic vortices is also addressed. In the absence of conclusive analytical results, vortex 
evolution is numerically investigated: the obtained results suggest that, contrarily to solitons 
with null angular momentum, vortex are unstable. 

1 Introduction 

Roughly speaking, a vortex is a solitary wave ip with non- vanishing angular momentum M (?/>). 
A solitary wave is a solution of a field equation whose energy is localized and which preserves 
this localization in time. The vortices in the Klein-Gordon equation (KG) are also considered 
in the Physics literature with the name of spinning Q-balls even if they do not exhibit spherical 
symmetry (see e.g. [53], [T3]L 

We recall also some existence results of solitary waves and vortices for KG: 

• For the case M (-0) = 0, we recall the pioneering paper of Rosen [20] and [14], [22], [10] . 
The spherically symmetric solitary waves have been called Q-balls by Coleman in |15| and 
this is the name used in the Physics literature. 

• Vortices for KG in two space dimensions have been investigated in [19]; later also three 
dimensional vortices for KG have been studied (see [9], Q], [23], [l"3].[6]L 

The aim of this paper is twofold. First, we give some existence results in situations not 
considered in the literature, namely we prove the existence of vortices in the case N > 3. For 
the N = 2 case see [8]. We introduce a new method, based on the energy minimization, which 
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allows to find vortices with a prescribed charge in a suitable class of functions (charge is defined 
in Section [20. 

Second, by taking advantage of this method, a numerical investigation is presented. In 
particular, some results in dimension N = 2 are reported, namely exemplificative vortex profiles 
by varying charge and angular momentum, together with relevant trends for vortex frequency 
and energy-charge ratio. Also, we discuss the stability of these vortices; such a problem is 
currently open, but we present some numerical simulations which suggest that the vortices 
arc unstable. In particular we show that stability cannot be proved by using the standard 
methods. In fact, in order to establish orbital stability, usually one uses the fact that solutions 
are minimizers of the energy functional on the manifold of fixed charge. Unfortunately it turns 
out that this is not the case, as we show in Theorem 12.121 actually our vortices are minimizers 
only in a suitable class of functions. 

This paper is an extension of [5], where we addressed the same problem for solitary waves 
with M (V>) = 0. 

Remark 1.1. In many situations, as in this paper, the existence of stable structures such 
as solitary waves and/or vortices is obtained by minimising the energy (which is defined in 
Section [2]) over a class of configurations of a given charge. If such a minimizing configuration 
exists, we may think that there is a force binding the "matter" (see [7] for details). The 
corresponding solitary waves have been called hylomorphic in [5] (see also [8]), by merging the 
Greek words "/ij//e" = "matter" and u morphe" = a form" . For this reason, the solutions considered 
in the present paper can be called "hylomorphic vortices" . 

The paper is organized as follows: in Section [2] wc prove the existence of vortices in the case 
N > 3 under general assumptions. Moreover the stability problem for the vortices solutions to 
(|NKGI) is addressed in subsection 12.21 Eventually the numerical investigation is presented in 
Section [3] 



2 The nonlinear Klein-Gordon equation (NKG) 

Let us consider the nonlinear Klein-Gordon equation (NKG) 

□V> + W'(tp) = (NKG) 
where □ = d? - A, ip : R N ->• C with N > 2, and W : C -> K with 

WW = Fm, W'^)=F'M)^ (2.1) 

for some smooth function F : M + — > ML Equation (|NKG[) is the Euler-Lagrange equation of the 
action functional S = J Cdxdt with Lagrangian density 

£(V>,£W,VV0 = \ |9 t V'| 2 - ^|VVf - WW). 
It is useful to write ip in polar form, namely 

iP(t,x) =u{t,x)e lS{t ' x) 

where u(t, x) <G R + and S(t,x) € M/(27rZ). If we set u t = d t u and St = d t S, the state ip is 
uniquely defined by (it, u t , 5 t , Vu, VS). Using these variables, the Lagrangian density C takes 
the form 

C{u,u t ,S t ,Vu,VS) = - u\ - |Vw| 2 + {S 2 t - |VSf) u 2 - F{u) 



1 While we were writing this paper, a similar existence result appeared in [3] for the N = 3 case under slightly less 
general assumptions. 
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and equation (|NKG[) is equivalent to the system 

□« + (|VS| 2 - S 2 ) u + F'{u) = 

d t (u 2 S t ) + V • (u 2 VS*) = 0. 

Nocthcr's theorem states that any invariance for a one-parameter group of the Lagrangian 
density implies the existence of an integral of motion, namely of a quantity on solutions which 
is preserved with time. Thus equation (|NKG[) has for N = 3 ten integrals: energy, momentum, 
angular momentum and ergocenter velocity. Moreover, another integral is given by the gauge 
invariance: charge. Easy computations show that the integrals of motion have the following 
expressions with respect to the variables (-0, dtip, VVO or (u, u tl St, Vu, VS): 



• Energy: 



£ =J i|9^| 2 + i|V^| 2 + ^(V0 



dx 



• Momentum: 



= / [| {u t f + \ |Vu| 2 + i [S? + |VS| 2 ] u 2 + F(u) 
3 = -Re J dtip\hp dx = - J (u t Vu + S t VS u 2 ) dx 



(2.2) 



dx 



• Angular momentum: 

M = Re J x x Vi&op dx = J (x x V S S t u 2 + x x Vu u t ) dx 



(2.3) 



Ergocenter velocity: 



K = tP 



= tf— Ix 



i|9^| 2 + i|V^| 2 + ^) 



dx 



Charge: 



H = Im / dtipip dx = St u 2 dx 



(2.4) 



A method to obtain vortices or spinning Q-balls is to write x £ K w as x — (y, z) € M 2 x M. N 2 
and to look for solutions of the form 



ip(t,x) = u{x) e l(i0{v) -^\ 



u>0, wei leZ 



(2.5) 



where 

%) = Im log ( Vl + iy 2 ) g M/2ttZ 

is the angular variable in the plane (2/1,2/2), for which we set r 2 := ?/ 2 + y 2 . With this ansatz, 
the NKG reduces to 

-Au+(£ 2 |V6»| 2 -o; 2 )u + F'(u) = (2.6) 
uA8 + 2Vu ■ V6» = 0. (2.7) 

By definition the function 9 satisfies 



A0 = 0, 
and if we assume that 



V0= - 



(-g,^,o,...,o), |V*|-i 

u(x) = u(r, z) 
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then (|2.7p is satisfied. For solutions of the form (|2.5|) . energy (|2.2[) and charge f|2 .4[) become 



E(u,u,e) = 



Jiv^l 2 



— 



u 2 +F(u) 



dx 



(2- 



H(u,u) = - J oju 2 dx (2.9) 
Moreover, in the particular case iV = 3, the angular momentum (|2.3[) becomes 

M(u,w,f) = (0,0,- / £wu 2 dx) = (0,0,£H(u,u)) (2.10) 



which is a non- vanishing vector if I ^ and H(u,lj) ^ 0. For this reason solutions (|2.5|) are 
called vortices or spinning Q-balls. 



2.1 Existence results for vortices for N > 3 

Letting x = (y, z) £ R 2 x R w_2 and r := yl/f+yf we consider the space of cylindrically 
symmetric functions of the form u(x) = u{r,z). In particular, letting O be an open subset of 
R^ -2 we consider the Hilbert space 



Hl(R 2 xO) = Mx) = <p(y,z) : <p £ C °°((R 2 \ {0}) x 0)} n ' nB ' 

obtained as the closure of cylindrically symmetric functions in C^°((R 2 \ {0}) x O) with respect 
to the norm 



IVtil 



1 



1 



dx. 



We denote by j| • \\fji and || • ||2 the standard norms of the spaces i/ 1 (R Ar ) and L 2 (M. N ). 

We remark that the cylindrical symmetry is sufficient to recover enough compactness thanks 
to the following lemma which is an immediate consequence of a result of Esteban-Lions |17j 

Lemma 2.1. [Compactess Lemma] Let O a bounded open subset ofM. N ~ 2 . 
H^(R 2 X0)4 L P (R N ) is compact for all p £ (2, 2*). 



Then the embedding 



As discussed above, we consider the cone of non- negative functions u{x) which are cylindri- 
cally symmetric and in H^(M. 2 x O) with O = 1 N_2 . Hence we introduce the notation 



Hi :={ueHl{i 



oN-2 



) : u(x) > j 



We are thus led to prove existence of solutions (u, lu, £) £ Hi x R x Z to equation ()2.6|) with Unite 
energy (I2.8[) . Since we are interested in spinning solutions, we fix I ^ 0, and look for solutions 
(it(£), u>(£)) with oj(1) ^ 0. In the following we drop the I dependence in notation. Without loss 
of generality we can restrict ourselves to the case u> < (so that charge is positive). 

First of all we remark that a variational approach to equation (|2.6j) in the space Hi x R x Z 
gives the existence of solutions in a weak sense. Then, following for example the argument 
of Theorem 2.3 in [B], one can prove that any weak solution is also a solution in the sense of 
distributions. Hence we can look for weak solutions to (|2.6p . Our variational approach is based 
on the following proposition, which is proved by a standard Lagrange multiplier argument (see 
[4] for an analogous result) 

Proposition 2.2. Let £ £ Z\ {0} be fixed. A couple (u,ui) £ Hi x R~ is a weak solution to 
V2. 6\) if and only if it is a critical point for the energy E(u,u>,t) constrained to the manifold 
Ch := {{u,uj) £ Hi x R~ : H(u,uj) = H(u,uj) = h ^ 0} of fixed charge. 
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We have thus obtained a very simple criterion to prove existence of vortices, we just need to 
prove the existence of constrained critical points for E to Ch- This problem can be solved by 
adapting existence results in [6] to this case, along the ideas introduced in [4], or by applying 
more general results from [8] . However we now give a novel proof, which permits to numerically 
construct the vortices as described in Section [3] 

Let us discuss the assumptions on the nonlinear term W in (|NKG[) . We write W as 

2 

777 

W(i/>) = F(\i>\) with F(s) = — s 2 + N(s), s>0 (2.11) 

and assume that 

• (W-i)(Positivity) F(s) > 0; 

• (W-ii)(Nondegeneracy) F(s) is C 2 with F(0) = F'(0) = 0, F"(0) = m 2 > 0; 

• (W-iii)(Hylomorphy) there exists Sq > such that N(so) < 0; 

• (W-iv) (Growth condition) at least one of the following holds: 

— (a) there are constants a, b > and 2 < p < q < 2N/ (N — 2) such that for any s > 

\N'(s)\ < as?- 1 +bs q -\ 

— (b) there exists s\ > sq such that N'(si) > 0. 

Remarks 2.3. We make some comments on assumptions (W-i), (W-ii), (W-iii), (W-iiii). 
(W-i) By (|2.8[) . (W-i) implies that the energy is positive. 

(W-ii) The necessary condition for the existence of solitary waves for (|NKG|) is F"(0) > 0. 
Results for the null-mass case, F"(0) = 0, are obtained in e.g. [10] and [2]. Here we consider 
only the positive-mass case. 

(W-iii) This is the crucial assumption which characterizes the potentials which might produce 
concentrated solutions. As we will see, this assumption permits to have states ip with hylomorphy 
ratio A (ip) < m (see f|2.13[) below). 

(W-iv) (a) This assumption contains the usual growth condition at infinity which guarantees 
the C 1 regularity of the functional. If we assume alternatively (W-iv)(b), the growth condition 
(W-iv) (a) can be recovered by using standard tricks (see [4]). 

Our main existence results is the following 

Theorem 2.4. Under assumptions (W-i)-(W-iv) and for any fixed I S Z \ {0}, there exists 
ho G K + such that for all h > ho the nonlinear Klein- Gordon equation (|NKG| admits a vortex 
solution ip(t,x) of the form (|2.5|) with finite energy (|2.2j) . charge (|2.4[) H = h and, for N = 3, 
non-vanishing angular momentum 



For fixed £ G Z\ {0}, it is sufficient by Proposition [221 to show that the energy E(u,u>,£) has 
a point of minimum on the manifold Ch for h big enough. However, it is possible and useful for 
numerical aims to reduce the problem to the minimization of a one- variable functional. Namely, 
for fixed I G Z \ {0}, for all the couples (u,w) G Ch the energy functional E{u,oj,1) can be 
rewritten as 



J h (u,£) :=E(u,Lo(u,h),£) = 



l\Vu\ 2 + ~u 2 + N(u) 



2 



where £ and h are parameters. The existence of a minimum for E on Ch is then equivalent to 
the existence of a minimum of Jh on \ {0}. 

Finally we introduce a fundamental tool in our variational approach, the quantity 

a, \ E(u,u>,£) 

A(u ' w): =i7(^r ( } 
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that we called the hylomorphy ratio (see [S]). 

Theorem 12.41 follows from the two following results 



Theorem 2.5. Under assumptions (W-i)-(W-iv) and for any fixed £ € Z\ {0} ; if there exists 
a couple (u,uj) € Cn such that A(u,cj) < m then Jh has a point of minimum on H\ \ {0}. 

Theorem 2.6. Under assumptions (W-i)-(W-iv) and for any fixed £ € 7L \ {0}, there exists 
ho 6 R + such that for all h > h^ 



inf A(w, uj) < m 



(2.14) 



The existence of hylomorphic vortices for all charges can be obtained by a stronger version 
of assumption (W-iii) . Let us consider the condition 

• (W-v) (Behaviour at s = 0) N(s) < -s 2+£ with < e < | for s e R + small enough 
(N(s) is defined in [2TTT|) . 

Theorem 2.7. Under assumptions (W-i)-(W-ii)-(W-iv)-(W-v) and for any fixed £ <E Z \ {0} 
and for all h € R + the nonlinear Klein-Gordon equation (|NKG|) admits a vortex solution ip{t,x) 
of the form (|2.5[) with finite energy (|2.2|) . charge (|2.4j) H = h and, for N = 3, non-vanishing 
angular momentum (|2.3[) . 



We first give the proof of Theorem 12. 5[ which is based on some preliminary results. 
Lemma 2.8 (A priori estimates). Assume conditions (W-i), (W-ii) and (W-iv)(a), and let 

m(£,h) := inf Jh{u,£). 

u£Hl\{0} 

Then there exist K\, Ki > such that if (u„) C H\ satisfies 

J h (u„,£) -> m(£,h) 
then K\ < ||w n ||2 < ||w n j|^fi < ||itn||^i < K 2 . 
Proof. First of all we can rewrite Jh(u,£) defined in (|2.12|) as 

"l ,_ i£ 2 



Mu,e) 



J 1 h* 

dx + — 



2 WuU 



(2.15) 



which is thanks to (W-i) a sum of non- negative terms. Hence any minimizing sequence (u„) 
satisfies 

f u 2 

\\Vu n \\ 2 < ai, /-jdx<a 2 , \\u n \\ 2 > a 3 

for some positive constants 0:1,02,03. 

In order to prove that ||M n ||^i < K2 it remains to show that || w„ || 2 < 04 for some positive 
CK4. Here we follow an argument from U which we include for completeness. From (W-ii) it 
follows that 

35 > 3ci > 0, such that F(s) > c x s 2 for < s < 5. (2.16) 
Let us assume by contradiction that 

\\unh °°- 

By ([2~T5| / F(u n ) dx is bounded and by (W-i) and ([2~T6| 



F(u n ) dx > 



0<u„<S 



F{u n ) dx > c\ 



i„ dx. 



0<u„<5 



(2.17) 
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On the other hand 



thus we have by (|2.17l) 



\\u n \\l 



I u 2 n dx + \ u 2 n dx — >• oo, 

Ja<u n <8 Ju„>5 

u 2 dx ^ oo. 



>>5 



This drives to a contradiction since for 2* 



1 



£2»-2 

and by the Sobolev embedding theorem 



2N 
JV-2 



it dx > 



u n >5 



ii dx 



u n >5 



u n dx < I U l n dx < C2 ||VM n || 2 2 < C2 «1 



,><5 



2* 



where C2 is the Sobolev embedding constant. 

Lemma 2.9. Let u n be a bounded sequence in H]. such that 



□ 



> 8 > for some q £ I 2, 



27V 
N - 2 



Then, up to subsequence, there exist a sequence (r n ) C 2 and it € u ^ 0, smc/i i/iai 
Un(" + r n) converges weakly to u. 



Proof. See Lemma 3.6 in [6]. 



□ 



Finally we need the following compactness result shown in [5] for a constrained minimization 
problem of functionals on the L 2 balls 

Lemma 2.10 (Minimization problem on L? constraint [6j). Let us consider the following min- 
imization problem 



G p := inf G{v) 



v£B„ 



where 



G(v) 



1£ 2 

2 1 ' " 1 ' 2 r 2 



\Vv\ 2 + -^v 2 + N{v) 



dx, 



and Bp = {u£ : \\v\\\ = p}. Under assumptions (W-ii), (W-iii) and (W-iv), if there exist 
vq £ Bp such that G(vq) < then for any sequence v n such that 

G(v n ) -> G p 

there exist, up to subsequence, a sequence (r ra ) C R w ~ 2 and av£ H\, v ^ 0, such that v n (- + T n ) 
converges to v in the H\ norm. 

Proof of Theorem \2.5[ First of all, fixed £ £ Z\ {0}, we can rewrite the hylomorphy ratio (|2 . 1 3[) 
as a function of the single variable u, that is 



J h {u,£) 1 



1£ 2 
2r 2 



\Vu\ 2 + -—u 2 + N(u) 



dx 



2 VlkHI h 



(2.18) 
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and remark that for all u £ H\ and all ft. > 

if h m 2 ,. Il9 \ . _ . 

2{Hl + - Ml )- m - (2 - 19) 

Therefore, if there exists (u,uj) £ Ch such that A(u,lu) < to, by (|2.18[) there exists u £ Hi such 
that Jh(u,£) < mh. Hence, thanks to (|2.19|) . for any minimizing sequence (u n ) C Hi such that 
Jh(u n , €) — ¥ m(£ 7 h) there exist n and eo < such that for all n > n 



G{u n ) 



1 I 2 

2 l Vu "| 2+ 2^ W " + iV(u " ) 



dx < £q < 0. 



This fact and the growth condition (W-iv) guarantee that there exist S > and q £ ^2, j 2 ^ 

such that 1 1 u n 1 1 > 8 for n > n. Hence we can apply Lemmas 12.81 and 12.91 to conclude that up 
to subsequence there exists a sequence (t„) C R*" 2 such that 

«„(■ + T n ) -^u ^0 weakly in Hi (2.20) 

with 

J u 2 n dx -> p > (2.21) 

by Lemma \2. 81 We claim that 

to 2 h 2 

m (e i h) = G p + —p+- (2.22) 

where p is defined in ()2.21j) and G p := mis G (see Lemma fe.lOj) . In order to prove (|2.22j) we 
take v n = ||^| 2 u n G Bp and show that 

G(v n ) - G(u n ) -> 0. (2.23) 

Indeed, 

|GK) - = - l)i /[|V U „| 2 + ^i 2 ]dx + jN( J ^ r2 u n ) - N(u n )dx\ < 

< KpfM - !)l5 /[|V«n| a + £« 2 ]^ + l(^ - 1)1 / \N'{6^-u n + (1 - 6)u n )\u n dx < 



< 



for some function 9 : 1 1 *' — > (0, 1). Moreover, ||un|| ^j-i is bounded by Lemma \2~El and applying 
(W-iv) (a) we get 

/ WiOj^Un + (1 - 9)u n )\u n dx < 



< a 



and the right-hand side is bounded thanks to the Sobolev embedding theorem. Hence (|2.23[) 
follows since ,, P M , — > 1. Hence it follows that 



%)4m(a)-yP- 



S 



and 

MG< m (£,h)-^p-f. 
B p 2 2p 



The opposite inequality holds since 



1 h 2 m 2 „ 2 \ .„,,„, m 2 h 2 



inf G(u)= inf [ J h ( u ^) ---— + — \\ u \\i = inf J h (u,l)-— p-— > 
ueB p ueB p \ 2 \\u\\2 2 J ueB p 2 2p 

.. , . ttl 2 h 2 
> m (l,h)-—p-- 

where in the last inequality we used the fact that m(£, h) is computed on a set containing B p . 

Hence we have proved that any minimizing sequence (it„) C H\ for J^. with m(£, h) < mh, 
satisfies (|2.20|) and (|2.21j) . and we find a sequence v n £ B p which is minimizing for a functional 
G satisfying the assumptions of Lemma 12.101 Indeed from (|2.22|) it follows that 



m 2 h 2 

-7TP+7T 
2 2p 



G P + — P+ — <mh 



whence that G p < 0. And from (|2.23|) it follows that G(v n ) — > G p . Hence, applying Lemma 
12.101 it follows that, up to subsequence, there exist a sequence (t„) C R n ~ 2 and a v e Hi, 
v ^ 0, such that v n (- + r n ) converges to v in the Hi norm. Hence 

u n (- + r n ) = ^""jj 2 !)„(■ + r„) ->• v in the Hi norm 

since ||u n ||2 \fp~- Theorem 12.51 is proved. □ 
Proof of Theorem \2.6l We start the proof by defining 

a := inf (2.24) 

sGR+ ks 2 



Thanks to (W-iii) we have 



ao < m 2 . 



The hylomorphy ratio (|2.13p can be written also as 

E(u,u) 1 1 
H(u,ui) 2 2ui 

where 

_ l(\\Vu\ 2 + \f\u\ 2 + F{u))dx 
J I it 2 dx 

Hence, for any fixed it, we have 



we 

Hence it is sufficient to show that 



inf A(u, lo) = A (u, -Ja(u)\ = \l a(u). 

'EM - V / 



inf a (it) < ao < m . 
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By definition of a>o (|2.24p . for any fixed e > there exists sq E ffi+ such that 



12 <"o+2 

2*0 Z 



Then let us consider the sequence of functions u n (x) = u n (\y\,z) = f{\z\) v n {\y\) with 



for \y\ <R n -l 

s (\y\ - Rn + 1) for Rn - 1 < \y\ < R n 



v n (\y\) = < 



so 



for R n < \y\ < 2R r , 



/(N) = 

and assume R n —> oo. Then u n £ Hi and 



s (2R n - \y\ + 1) for 2R n < \y\ < 2R n + 1 
for \y\ >2R n + l 

1 for < \z\ < 1 

2-\z\ for 1 < \z\ < 2 

for \z\ > 2 



-u n dx 



2-R„ «2 



J \Vu n \ 2 dx = 0(R n ) 

f(z)dz) + (ln(i?„)) = OQn(i? n )) 



-rsl dr 



2R„ 



\z\<2 



\ I rF(sa)dr) (J f(z)dzj+o(R n )< 



< 



(C 


1 

r — 




2 


j\\Un 


2 dx 



' 2' 

2R„ 



f(z)dz + o(R n ) = const(a Q + -)s R n + o(R n ) 



\z\<2 



r-sldr 



f(z)dz + o(R n ) = const s R n + o(R n ) 



\z\<2 



It follows that a(u n ) < a$ + e for sufficiently large n. Moreover 

H (u n , -\fa(u n )^j « R n . 

We have thus proved that there exists ho £ K + such that 

inf A(u, tu) < m. 

Let now h > ho and (u,lj) £ Ch a such that A(u,lj) < m. We want to show that there exists 
(u,ui) <G Ch such that A(u,w) < m. Recalling that x = (y,z) 6 R 2 x M w_2 , let us define for 

u\(y,z) = u j ^ = w 

Then (ua,w) € Ch if 

— cd \ \u\W 2 . = —qj A 2 II Oil 2 = — uj A 2 — ^- = h 

—lu 
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that is, if A 2 = j~ > 1. Moreover, for the single terms appearing in A it holds 

j\Vu x \ 2 dx = J (|V y w A | 2 + |V 2 w A | 2 ) dx = J {\V y u\ 2 + \~ 2 \V z u\ 2 ) dx < J \Vu\ 2 dx, 



since r 2 = \y\ 2 , and 
Hence using also ft. > ho it follows 



— u x dx = I — it dx 



N{u x ) dx = A 2 / N{u) dx 



k(u\,w) 



< 



1 

~ h 
1 f 1 



|v^| 2 + -- Wa + tvK) 

2 r z 



dx 



1 

fro 



1 1 



A 2 



2 V Kill ft 



"A 9 < 



dx + — N(u) dx + - [-co 



2r 



and the proof is finished. 



□ 



Proof of Theorem \2.7\ Arguing as in the proof of Theorem 12.41 we only need to show that for 
any fixed £eZ\{0} and for all ft e R+ 



inf A (it, uj) < m 



(2.25) 



Choosing uj = —m we have 



inf A(u,u>) < inf A(it, — m). 



Notice that A(u, — m) < m if 



I| Vu | 2 + i^| M | 2 + iV( u ) 



Hence it is sufficient to show that for any p > 



inf 



\\^ u \ 2 + 5^2 H 2 + N{u) 



where 



Hl p :={ueHl s.t |M| 2 = p}. 



Now let us consider the sequence of functions u n {x) = u n Qy\,z) = f(\z\) v n (\y\) with 



for \y\ < R 
j£\y\ for R„ < \y\ < 2R n 



M\y\) = < 





s n for 2R n < \y\ < 4i?„ 

5s»-j£\y\ for 4i?„ < \y\ < 5R„ 

for \y\ > 5i?„ 
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1 for < \z\ < R n 

l-^-d^l-i^) for R n < \z\ < 2\R n 
for \z\ > 2Rn 

and assume R n — > oo and s n — > such that J \u n \ 2 dx = p. Notice that 

\u n \ 2 dx=( [ v n (\y\) 2 dy] ( [ f(\z\) 2 d 



such that a simple scaling analysis shows that lim„^oo s^i?„ = 7 > 0. Moreover we get 
[ \Vu n \ 2 dx=([ \Vv n \ 2 dy)(f f(\z\) 2 dz) + ( f v n (\y\) 2 dy) ( f \Vf\ 2 dz)=0(R 

JR N \JR 2 / \JR N - 2 / \JR 2 / \JR N - 2 J 

f N{u n )dx < - f \u n \ 2+t dx = -( f v n (\y\) 2+e dy) ( f f(\z\)*+*dz) = -0(R~ ieN ) 

JR« JR<* \JR 2 J VJR"- 2 / 

/ ^rjXJ r v n(\v\fdyj (j[ N 2 f(\z\) 2 dz^ = o(i?,- 4 ) 

The previous estimates imply for R n 00 and < e < that 

'1- .2 , 1^, ,2 



that 



2 ivw„r + 2^i u "i + ^f"") ) dx °" 



inf A(m,w) < inf A(u, — m) < to. 



□ 



2.2 The stability problem 

We now discuss the possibility that the vortices solutions to (jNKGp found in the previous section 
are orbitally stable, in the sense of the following definition (see e.g. |18j and 0]). For simplicity 
of notations we consider the case N = 3, so that we can use expressions (|2.3p and (|2.10j) for the 
angular momentum. However the same argument holds in any dimension N > 3. 

Definition 2.11. A vortex solution ^(t, x) = u(x)e< M{ ^- ut ^ is called orbitally stable if the set 



>, £, w) := {u(x + x ) e '<*%+©—*-») : x = (£, C ) e 



£ 2 x 



is stable in the -ff 1 x L 2 norm, i.e. for any e > there exists 5 > such that d(ip(0, < S 
implies d(ip(t, ■), T) < e for all £ € R, where 

d(^(t, ■), •)) := ||VC*» ') - u{-)e< l ^~^\\ m + \\8Mt, ■) - (-icj)u(-)e< M U-"Q\\ L z (2.26) 

To establish orbital stability there are essentially two methods used in literature, one based 
on the Lions's Compactness-Concentration Lemma and one introduced by Shatah (see [18j). 
However both methods use the idea that stability can be obtained for solutions that are non- 
degenerate points of global minimum of the energy functional constrained on the manifold of 
fixed charge. Hence we could investigate whether our solutions $>(t,x) = uip^e 1 ^^ - ^ are 
points of minimum for the energy £ constrained on the manifold of fixed charge %(V>) and fixed 
angular momentum M(?/>). Unfortunately it turns out that this is not the case, as we show 
in the main result of this section Theorem 12.121 However it remains open the possibility that 
these solutions are non-degenerate points of local minimum, as it happens for solutions found 
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in for (|NKG[) for solitons with vanishing angular momentum. This weaker property would 
be sufficient for orbital stability, but this is still an open question. 

Let us recall the expressions £ (|2.2[) and E l|2.8p for the energy on a general function ip £ 
H 1 (M 3 ,C) and on the vortices solutions satisfying the ansatz (|2.5j) respectively. The solutions 
ip in the previous section have been found as points of global minimum for the energy E on the 
manifold of couples (u, uS) with fixed charge. We now show that these solutions do not attain 
the global minimum of the energy £ . 

Theorem 2.12. For any fixed i £ Z\ {0}, letip(t, x) = u(x)e*( w M _£M ) be a solution of (|NKG[) . 
with h — T-L(ip) and m = M(^), found as in the proof of Theorem \2.4[ that is E(u, Q) is the 
minimum of E on the manifold Ch ■ Then 

inf £ = inf £ = inf E(u, u, 0) < 8($) 

H=h,M=m U=h {u,u)&C h 

Proof. Step 1. £($) > inf E(u,u,0) 

{u,u:)eCh 

First of all it holds 

£(tp) = E{u, u, i) > inf E(u, u, I) 
(«,w)eCa 

which is given by construction. Moreover we show that 

e r .i := inf E(u,lu,£) > inf E(u,ui,0) =: e 
(u,ui)ec h (u,u))ec h 

Let («„,«„) £ Ch, with u n £ H^, be a minimising sequence for e r ^. By Lemma 12.81 there exist 
Kx,K 2 > such that 

K x < \\u n \\ L 2 < \\u n \\ H i < \\ujjj, < K 2 (2.27) 
If by contradiction e r ^ = e then, using the notation -Br(O) := {y £ R 2 : r < i?}, we have 

lim L 3 ^dx = Q (2.28) 

lim f u 2 cfa; = Vi? > (2.29) 

Indeed by definition of E(u,w,i) it immediately follows that 

1 f I 

e r f >e+ inf - / —u 2 dx 

(u,L0)£C h 2 J R3 H 

hence (|2.28|) is proved. Let us now write for any R > 

u 2 dx J = 4tt 2 ( / ( / ru i n dr\dz\ < 

<4^ 2 [ / [ R r 3 u 2 n drdz) ( [ 1^ ^ drdz] < An 2 R 2 \\u n \\ 2 L2 [ % dx 
\JrJo J \JrJo r J Jrz r 

Hence (|2~29|) follows from (|2~27f and ([2^28)1 . By ((2T29|) it follows that u„ converges to weakly 
in L 2 , which contradicts (|2.20l) (see also Lemma [ 

Step 2. inf £ = inf E(u,w,0) 

H=h (u,u)eC h 

See g]. 

Step 3. inf £ = inf £ 

« = h.M=m H=h 



\ \ 2 

,2 
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The inequality 

inf £ > inf £ 

U=h,m=m U=h 

is immediate since on the left hand side the infimum is taken on a smaller set, since we constrain 
also on the manifold of fixed absolute value for the angular momentum. The opposite inequality 
is proved by constructing a sequence {4>n} such that H(ip n ) = h and M(^>„) = m, which satisfies 

£Ub n )^r inf E(u,u,0) = inf £ (2.30) 

For any fixed e > let U(x) and v(x) be non- negative radially symmetric functions in (R. 3 , R + ) , 
let ui < and £ £ Z which satisfy 

U{U{x)e- wt ) = - J RN ujU 2 dx = h-e (2.31) 

E{U,u),Q)= inf E{u,u,Q)+E (2.32) 
(ti,w)ec h _ e 

M(U(x)e~ %ut ) = (2.33) 

U{v(x)e< 19 ^-^) = - J RN ujv 2 dx = e (2.34) 

M(v(x)e« a M- u *)) = (0, 0, -£ J RN ujv 2 dx) = m (2.36) 

Recalling the notation x = (y, z) G M 2 x K, we now define 

^(s) = C/(x)e-^ t + v(y - y n , z ) e *Mv)-»*) 

where {y n } C R 2 is a sequence satisfying \y n \ — > oo. For n big enough the supports of U 
and v(y — y n , z) are disjoint, hence since the integrals in the expressions of charge and angular 
momentum are translation invariant, we obtain 

U{i> n ) =H(U(x)e- wt ) + H(v(x)e< ee M-^) = h 

from QOP and (12341 . 

M(V>„) = M.{U{x)e' lut ) + M(v(x)e< M ^- u ^) = m 

from (|2.33[) and (|2.36[) . We now compute the energy of tp n . Again, since supports are disjoint, 
we can write 

inf E(u, u, 0) < e(if> n ) = £(U(x)e-* ut ) + £(v(x)e z( - ee ^-^) = 
(u,u>)ec h 



1 



E(U,u,0)+E(v,w,0) + ^- I ^v 2 (y-y n ,z)dydz 

2 7b3 r 2 



r- 



2 



if £ 

inf E(u,uj,0) + e + E(v,uj,0) + - / v 2 {y,z)dydz< 
(u,u)ec h - e 2 7 R3 (r + |y„|)^ 

< inf E(u,uj,0) + e + is(2j—)+- — \\v\\ 2 L2 

where z/(-) is the modulus of continuity of the energy E{u, ui,Q) with respect to the H 1 norm, 
and we have used (|2.32[) and (|2.35[) . 



14 



Finally we remark that 



inf E(u,uj,0)= inf 
(«,u)ec k - e «Gff 1 \{o} 



l -\Vu\ 2 + W(u) 



< inf 



i|Vu| 2 + VF( U ) 



1 h 2 



2 J u 2 dx / (u,w)eCh 



l (h-er 



inf E(u,u,0) 



We can conclude that there exists n(e) for which 

inf £(u,w,0) <£(t/w<o) < inf w, 0) + 2e + v { 2. 

(u,to)eCh (u,ai)eC ft \ ' 

and (|2.30p is proved. 



□ 



3 Numerical approach to hylomorphic vortices for NKG 

In the current section we consider for ease of presentation only the case N = 2, so that u(x) = 
u(r) is a radially symmetric function. Existence of hylomorphic vortices for N = 2 is shown in 
[8]. In particular, once introduced a straightforward method for the numerical construction of 
vortices, we present relevant simulations addressing their orbital stability. 



3.1 Numerical construction of vortices 



It has been shown in Section [2TT1 that, for fixed I € Z\ {0}, vortices can be found as points 
of minimum for the energy E(u,ui,£) on the manifold Ch of fixed charge H(u,uS) = h. Such 
a constrained minimization problem in two variables can be reformulated as an unconstrained 
minimization problem in one variable, as shown above in (|2.12j) . In particular, hereafter we can 
use u £ Hi as independent variable and, letting 



co (u, h) 



h 



J u 2 dx ' 

numerically study the minimization problem for the hylomorphy ratio 



Jh{u,t) 



\\^u\ 2 



1 £ 2 

--u 2 + W(u) 



dx 



2ju 2 dx' 



(3.1) 



(3.2) 



on H\. To this purpose, we consider the evolutionary problem (which generalizes the one treated 
in[l]) 



d T u(r, t) = -h dk h = Ait - W'{u) 
u{r,r) = 



I 2 



u(r, t) 







in (0,f) x K+ 

on {r = 0} x M+ 
on {r — f} x R + 



(3.3) 



in which r represents a pseudo-time, r denotes a chosen upper bound for the r— domain (dis- 
cussed below) and u> = lj(u, h) as in (|3.ip . The evolution of u according to (|3.3[) is a gradient 
flow and therefore a non-increasing trend of A/j(u(r, r)) versus r is obtained, well suited as the 
sought minimization process. 

The problem (|3.3j) was then discrctized by a classical line method; in particular, 2nd order 
and 1st order accurate finite differences were respectively used for space and time discretization. 
Moreover, f was chosen large enough to suitably contain the vortex profile support, and the 
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chosen charge h was directly enforced at each time level by evaluating the frequency w™ (super- 
script n = 0,1,2,... hereafter denoting the n-th time level) from the corresponding numerical 
solution through the discrete counterpart of (|3.1[) . Finally, time- advancing was stopped when 
\uj n+1 — u n \/uj n < e u , where represents a predefined threshold (a relative error on A/j was 
considered as well). It is worth remarking that: (i) the proposed method managed to efficiently 
converge by starting from several initial guesses; (ii) for many of the carried out numerical 
experiments it was necessary to adopt a very severe threshold, e.g. e u = 10~ u , in order to 
obtain discretization-independent results (this point suggested the possible presence of rather 
flat regions around the sought minimum, somehow challenging the achievement of the "true" 
radial profile). 

Let us now consider, for the sake of illustration, the following potential: 

W (V) = M - log (1 + |V|) = \ IVI 2 - \ M 3 + °(H 3 ), (3-4) 



for which NKG specializes to 



Figure [T] shows some radial profiles u(r) of two-dimensional (2D) vortices obtained by minimizing 
the functional (|3.2p coupled with (|3.4[) . through the procedure described above. It should be 
noticed that: (i) for a given £, the peak of the u profile significantly decreases by decreasing the 
charge h, while the associated support is less affected; (ii) for a given charge h, the peak of the 
u profile smoothly decreases and drifts away from the origin by increasing I. Moreover, Figure 
[3 shows characteristic trends for both frequency and hylomorphy ratio, still versus h and t. It 
should be noticed how both |w| and A tend to 1 for h —> 0, as well as for increasing values of I. 
Finally, once obtained u(r), it is straightforward to build vortices through (|2.5[) ; two examples 
are shown in Figure [3] (a grid is added for ease of rendering). 



3.2 Numerical experiments on the stability of vortices 

Numerical approximation of NKG has long been studied (see e.g. [23], [TH] and [IS]), in order to 
also investigate nonlinear phenomena like the dynamics/interaction of solitary waves or other 
coherent structures. In light of the relevant literature, we decided to preliminarily assess the 
suitability of a line method approach based on centered finite-differencing for Laplacian dis- 
cretization and leap-frog time-advancing, by virtue of its relatively simple coding. In particular, 
we firstly considered the time-evolution of a (non-rotating) hylomorphic soliton for which orbital 
stability is proved in [5], on a square domain with periodic boundary conditions (2D-torus). The 
adopted numerical scheme is symplectic, thus being able to accurately capture the considered 
dynamics over a long time- interval (see e.g. [12], [16] and [21]). 

More in detail, we firstly chose h = 500 and we obtained the radial profile and the frequency 
of a non-rotating soliton as in [5J (i.e. by a minimization strategy like Q3.3p . yet with I = 0), with 
W defined as in (|3.4[) . We then introduced a square grid with spacing 8x (and side length large 
enough to properly contain the soliton support), and we endowed the aforementioned soliton 
with speed v = 0.5 along a grid axis, by Lorentz transform. Let us denote by {ip th (x, t) ,ipl h (x , t)) 
the obtained (translating) theoretical soliton. Furthermore, once introduced the time-step St, 
we defined the initial conditions as ip° := rp th (x, 0) and ip^ 1 ^ 2 := ^p t t h {x 1 —6t/2)), where ~ stands 
for sampling on the computational grid, and the 5t/2 time-shift accounts for time-staggering 
of the leap-frog scheme [21] , Therefore, after n time-steps the adopted scheme provides ib n w 
$ th (x,n5t) and Vt™~ 1/2 « ^\ h {x, (n - 1/2) St). Time-stag gering clearly affects the numerical 
approximation of those entities which simultaneously involve ip and i/j tl like e.g. all NKG first 
integrals; however, corresponding error can be kept contained by adopting small time-steps, as 
required by stability constraints (see below). In light of this aspect and with the main aim 



16 




17 




3 10 20 30 40 50 °" ?5 10 20 30 40 50 



(c) (d) 

Figure 2: Frequency (a,c) and hylomorphy ratio (b,d) of 2D vortices: as a function of h, having 
fixed 1 = 1 (a,b); as a function of I, having fixed h = 700 (c,d). Circles in (a,b) are associated with 
the following values of h: {5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700}; circles 
in (c,d) are associated with the following values of £: {1,2,3,4,5,10,15,30,50}. W defined as in 
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(a) (b) 

Figure 3: Surface plot of Re(ip) for a 2D vortex with: h = 700, £ = 1 (a); h = 700, i = 4 (b) 
(corresponding radial profiles are shown in Figure [1]); W defined as in (|3.4p . 



of investigating orbital stability, a "staggered" discrete counterpart of (|2 . 26[) was introduced, 
namely the following orbital stability norm: 

M"-4> th (x,n6t)\\ H i , m- 1/2 -tt h (x,(n- 1/2) 6t)\\ L . 

o n q := = =— , n = 0, 1, . . . , (3.5) 

\\^(x,0)\\ H i M h (x,-dt/2)\\ L2 

where discrete H 1 and L 2 norms are tacitly understood (they were computed by classical fourth- 
order Simpson quadrature), and the denominators (involving time-conserved entities) simply 
act as scaling factors. It should be noticed that the adopted initial conditions enforce 5q S = 
0; subsequent evolution represents the main asset of the considered numerical investigations. 
Moreover, it is worth remarking that numerical approximations of NKG first integrals were 
defined similarly to ()3.5|) . yet corresponding expressions are here omitted for brevity. Formally, 
the considered numerical scheme introduces a truncation error 0(Sx 2 + St 2 ). However, stability 
constraints affect the adopted explicit time-advancing (see e.g. [24] and [21] for some model 
problems), so that a ratio St/Sx = 1/10 was adopted for the numerical experiments at hand 
(it was also checked that such a ratio provided time-discretization- independent results). Hence, 
it was space discretization to directly modulate the leading term 0(8x 2 ) of truncation error. 
Consequently, grid coarsening/refinement was regarded to as a mean for implicitly tuning a 
perturbation on the theoretical soliton profile which was exploited to seed the initial conditions. 
More in detail, several numerical simulations were carried out by varying Sx; corresponding Sq S 
trends are reported in Figure [4] 

Conservation of NKG first integrals was systematically monitored for all the considered simula- 
tions. For instance, energy variation relative to its time-average was in the range (—8, 8) • 10~ 2 , 
(—1, 1) ■ 10~ 3 and (—2,2) ■ 10 -6 , when respectively adopting Sx = 5, 1 and 0.1. Corresponding 
relative variations for charge were below 10 -12 in absolute value. These results strongly support 
the remarkable features of the adopted symplectic integration scheme, while also establishing 
sound foundations for the interpretation of the obtained Sq S trends. On regard, a rough dis- 
cretization like e.g. the one associated with Sx = 5 introduced a disrupting perturbation on 
the seeded profile: the numerical solution soon lagged behind the moving theoretical one. This, 
in turn, only periodically allowed for support overlapping: for almost all times, numerical and 
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2 



1.5 



1 

0.5 

°0 5 10 15 20 25 30 35 40 
t"/T 

Figure 4: Exemplificative trends of 5q S versus time t n = nSt for an orbitally-stable hylomorphic 
soliton (t = 0, non-rotating) with h = 500; W denned as in (|3.4|) . Time non-dimensionalised 
through the soliton period T. Each curve is associated with a grid; corresponding characteristic 
size 8x was varied in {5,2,1,0.5,0.2,0.1}. Non-rotating solitons were preliminarily simulated, in 
order to assess the suitability of the proposed numerical approach. 

theoretical solutions occupied disjoint regions of the computational grid, and this justifies the 
asymptotic value around y/2 in Figure |U Moreover, as grid spacing was reduced, both the 
initial solution was better sampled (with respect to the theoretical shape) and the truncation 
error decreased. This, in turn, incrementally led to grid-independent numerical results (see the 
curves clustering in Figure 21 where lower values of Sx are not shown for ease of readability), 
which also exhibit an orbitally stable character (the weakly increasing trend is merely due to a 
cumulative effect of discretization errors over time). The obtained numerical results, which are 
fully consistent with the underlying theoretical framework, encouraged to extend to the vortex 
case the proposed approach, as well as the idea of interpreting space discretization as a mean 
for issuing tailored perturbations with respect to theoretical profile of the considered solitons. 

The same numerical approach was therefore applied to rotating solitons, by firstly considering 
a vortex with £ = 4 and charge h = 500. In particular, a polar grid was introduced and the 
characteristic discretization size Sr along the radial direction was chosen as the leading term for 
truncation error. This was achieved by systematically adopting a finer discretization along the 
circumphercntial direction, and by choosing St/Sr = 1/50 in order to ensure stability (indeed, 
the ratio 1/10 previously adopted for non-rotating solitons turned out to be insufficient in the 
present case) . The initial conditions were defined as above, and the simulation was advanced by 
monitoring first integrals and the orbital stability norm defined by Q3.5p . In particular, several 
numerical experiments were carried out by varying Sr; corresponding Sq S trends are reported 
in Figure [5] (b) . 

Conservation of NKG first integrals was systematically monitored for all the simulated vortices. 
For instance, energy variation relative to its time-averaged value was in the range (— 1, 8) • 10~ 2 , 
(0,5) • 10~ 3 and ( — 1, 1) • 10~ 5 , when respectively adopting Sr = 5, 1 and 0.1. Corresponding 
ranges for charge were (—12, 1) • 10~ 2 , (—1, 1) • 10~ 6 and (—2, 2) • 10~ 10 ; corresponding ranges for 
the angular momentum were (—15, 5)-10~ 2 , (—1, 1) • 10 — 6 and (—2, 2)40~ 8 . These results further 
supported the adoption of the above described numerical scheme, and gave us confidence in the 
obtained Sq S trends. On regard, no stable dynamics was observed: vortex rupture occurred 
for all the considered test-cases, at the latest nearly t/T = 63 on the most refined grids (such 
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dr = 5 

-»-dr = 2 

2.5- ... dr=1 - 2.5- 

-*-dr = 0.5 
" '0 " dr = 0-1 
-a-dr = 0.05 

p _ ■ <"dr = 0-01 . p _ 




t n AT t n /T 

(a) (b) 

Figure 5: Exemplificative trends of 5q S versus time t n = n5t for hylomorphic vortices with £ = 4, 
having fixed h = 100 (a) and h = 500 (b); W defined as in (I3.4p . Time non-dimensionalised through 
the soliton period T. Each curve is associated with a grid; corresponding characteristic size 6r was 
varied in {5,2,1,0.5,0.1,0.05,0.01}. The obtained trends are remarkably different from the ones 
in Figure [H 



dr = 5 

-»-dr = 2 
.».dr=1 
-«-dr.0.5 

uO"dr = 0.1 
-•-dr. 0.05 
<"dr = 0.01 



a threshold is clearly visible in the figure). A rough discretization like the one associated with 
Sr = 5, for instance, immediately caused a phase lag between numerical solution and theoretical 
vortex, which in turn produced the corresponding periodic Sq S trend in Figured] (b). Moreover, 
perturbation reduction (as induced by grid refinement) did not prevent vortex rupture; most 
importantly, not even it affected rupture time (e.g. by introducing incremental delays). This 
aspect is most evident in the figure: rupture time was not significantly affected by a reduction 
of two orders of magnitude in the grid characteristic size (say, from 8r = 1 down to 0.01). 
Such results are completely different from the ones observed for hylomorphic solitons, i.e. in 
the orbitally stable case. Nonetheless, the aforementioned numerical experiments were repeated 
by considering a vortex with £ = 4 and charge h = 100, in order to possibly corroborate the 
preliminarily obtained results. Corresponding trends for 8q S are reported in Figure [5] (a): they 
fully confirm those observed for h = 500. This point, together with the fact that rupture 
occurred nearly t/T = 125 for such a smaller charge value, incidentally suggests that no trivial 
correlation between charge and orbital stability might be inferred at the present research stage. 
Overall, the carried out numerical simulations suggest that vortices may be unstable. 
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